Tailored particle current in an optical lattice by a weak time-symmetric harmonic 
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Quantum ratchets exhibit asymptotic currents when driven by a time-periodic potential of zero 
mean if the proper spatio-temporal symmetries are broken. There has been recent debate on whether 
directed currents may arise for potentials which do not break these symmetries. We show here that, 
in the presence of degeneracies in the quasienergy spectrum, long-lasting directed currents can be 
induced, even if the time reversal symmetry is not broken. Our model can be realized with ultracold 
atoms in optical lattices in the tight-binding regime, and we show that the time scale of the average 
current can be controlled by extremely weak fields. 
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Brownian motors or ratchets are spatially periodic sys- 
tems with noise and/or dissipation in which a directed 
current of particles can emerge from an unbiased zero- 
mean external force [l], Q • Models for biological engines 
that transform chemical energy into unidirectional me- 
chanical motion behave as Brownian motors Q. Ex- 
tensive studies of the ratchet effect in classical systems 
[1] stated the relation between symmetry breaking 
tcntials and the existence of the asymptotic current 
For a system driven by a flashing potential of the form 
V(x, t) = V x {x)V t (t) with V t {t) =V t (t + T) of zero mean 
and V x {x) = V x (x + L x ) there are up to four different 
symmetries in the classical system that must be broken 
in order to generate an asymptotic current [(| . A ratchet 
current arises if one breaks the relevant spatio-temporal 
symmetries, here denoted by S x , and the time-reversal 
symmetry St : (x,p,t) —> (x, —p, — t + 2t s ). Lately there 
has been an increasing interest in the coherent ratchet 
effect in Hamiltonian quantum systems 0- It has been 
shown that the same symmetry requirements apply to 
them Q, i.e. if the Hamiltonian preserves any of the 
symmetries, no asymptotic current is possible. 

Experimentally, directed current generation was first 
studied in solid state devices, quantum dots and Joseph- 
son junctions Q. More recently, the precise control 
achievable in cold atom experiments opened up the pos- 
sibility of realizing directed atomic currents for Hamil- 
tonian systems with controllable or no dissipation in the 
time scale of the measurements (HH. Recently, a very 
clean realization of a coherent quantum ratchet was ex- 
perimentally demonstrated in a Bose-Einstein conden- 
sate exposed to a sawtooth potential realized with an 
optical lattice- which was periodically modulated in time 
[13|. Directed transport of atoms was observed when 
the driving lattice potential broke the spatio-temporal 
symmetries. The current oscillations and the dependence 
of the current on the initial time and the resonant fre- 



quencies |14[ were measured, demonstrating the quantum 
character of the ratchet. 

Although the generation of an asymptotic directed cur- 



rent needs the breaking of the symmetries S x and St si- 
multaneously for unbiased potentials, there has been a re- 
cent discussion on the possibility of obtaining long-lasting 
directed currents without it 15-1|1]. Many-body effects 

15, with the proper choice of the initial state [H| 
or an accidental degeneracy in the quasienergy spectrum 
[l8| may result in a directed current without breaking the 
time-reversal symmetry. In contrast to previous works 
we show here that one can exploit a quasi-degencracy, 
present for a wide range of parameters, in the quasienergy 
spectrum in order to generate a long-lasting directed av- 
erage current in a weakly driven system where we can 
achieve full control over its magnitude and time scale. 
Previous work on quantum accelerators 12| has shown 
that in the presence of quantum resonances one can ob- 
tain large currents without breaking the time-reversal 
symmetry using a delta-kicked potential in time. Essen- 
tially, for particular values of the Hamiltonian parame- 
ters the spectrum of the Floquet operator becomes con- 
tinuous due to quantum resonances. Under such circum- 
stances one can obtain a linear increase in momentum 
with time which has been claimed as a true ratchet effect 
driven by resonances instead of noise [20j . However, there 
is no formal proof that the dynamics show unbounded ac- 
celeration for times longer that those that are computed 

12|. Nonetheless, a significant difference between quan- 
tum accelerator ratchets and our system is the existence 
of a constant component in the delta potential. 

One useful way of treating time-periodic quantum 
Hamiltonians, H(t) = H(t + T), is the Floquet formal- 



ism [22]. The cyclic states \<j>j{t + T)) = e 



are the eigenstatcs of the evolution operator for one 
period while the quasienergies £j are the eigenvalues. 
The solution to the time-dependent Schrddinger equa- 
tion H{t) | ip(t)) = ihd | ip(~b)) /9t can be spanned in the 
cyclic eigenbasis (H = 1) 



(i) 



where Cj = ((f>j(0)\ip(Q)). The average current generated 



2 



during n cycles is given by I(t = nT) = i ^m=i with 

i-mT 



X - I 



(ro-l)T 



(V»(<)|p|^(t))d*, 



(2) 



where p is the momentum operator. Note that due to 
the periodicity of the cyclic states the average current 
during n cycles can be simplified in terms of integrals of 
the cyclic states during the first period 



X{nT) 



(3) 



which is valid for 
where ip)jj' '■= 



a discrete Floquet spectrum and 

Our model, that represents well optical lattices, has by 
construction a pure point spectrum. In the limit of an 
infinite number of states the single-band tight-binding 
model remains integrable [231 ]. However, if more bands 
are added the spectrum could become absolutely contin- 
uous (in resonance) or singular continuous. Our guess, 
based on our previous experience in related problems [24| , 
is that except at very long times even if the spectrum is 
singular continuous the time evolution will be similar to 
that given by a pure point spectrum. In general, a sum 
of oscillatory off-diagonal terms with arbitrary exponents 
decays rapidly [25| and for long times only the diagonal 
terms in Eq. ([3]) remain. If both S x and St are broken, the 
cyclic eigenstates desymmetrise and carry net momen- 
tum, i.e. (p)jj ^0@. In such case, the asymptotic av- 



j.i ■ 



erage current at n — > oo is nonzero I(oo) = £^ ■ 
Correspondingly, if either of the relevant symmetries is 
not broken (p)jj = and thus I(oo) = 0. Note, how- 
ever, that the off-diagonal terms in Eq. ([3]) become rele- 
vant if the initial state projects mainly into degenerate or 
quasidcgcncratc cyclic states with (p) jj' 7^ 0. If one in- 
duces a resonance between the proper quasienergy states 
at low driving, the average current contains only a small 
number of terms in the sum and the exponents can be 
very small, leading to very slow oscillations whose pe- 
riod can be fitted by tunning the driving. In order to 
maximise the average current one should then optimise 
both the projection into the initial states 



■3^3' 



and the 



(p) jj> . We illustrate this here and show that it is possible 
to populate a high average momentum superposition of 
cyclic eigenstates for times which can be tuned up to the 
lifetime of the experiment. 

We consider a driven system of non-interacting bosons 
with H{t) = Ho + V(t) in a lattice of L sites with periodic 
boundary conditions |26[ and 



H = -J^|Z)(Z + 1| + |Z + 1><Z| 



(4) 



V(l,t) = Fsin(wt) 
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FIG. 1: (color online) a) Thick solid line shows the aver- 
age current, I(t) in Eq. ([3]), in recoil units (fcrecoii = 2ir/L) 
for an initial zero momentum state. Thin solid line is the 
current average per cycle T m in Eq. ((2j) obtained with a nu- 
merical integration of the Schrodinger equation for M — 5, 
L = 41, J/u = 1.0398, V/u = 0.1, (j>/iv = 0.2579 and a = 1.2. 
The driving frequency u is tuned to the resonance condition 
in Eq.©. Dashed line corresponds to the current of the ef- 
fective 3- level system in Eq. ([7]) and Eq. ((§]). The average 
current attains a maximum M. b) Numerical calculation of 
T(t) for same parameters except for solid line V/u — 0.1, 
= 0.2579; dashed line V/u = lO" 3 ^, 4>/k = -0.2579; 
dot-dashed line V/u = 10 -4 , cp/rv = 0.2579. Note that the 
time of peak average current scales with V~ 2 uj. 



with M integer, where J is the tunneling probability and 
1 1) represents the state of a boson located on site The 
eigenenergies of Hq are Ek = — 2 Jcos(27rfc/L) with inte- 
ger k = [-fc max , fcmax] where fc max = (L — l)/2 for L odd 
and the corresponding momentum eigenvectors (l\k) = 
e b ' 2 c L /y/L are degenerate for ±fc. For convenience one 
can introduce the basis (I \ Sk) = \/2/L cos(fc27rZ /L) and 
(I | dfc) = \JljL sin(fc27r//i) which are symmetric or an- 
tisymmetric under the inversion of k. We add a time 
and spatially modulated periodic function V(l,t) with 
frequency w = 2tt/T tuned to the M-dependent resonant 
condition El H3 



2lo = E 2 M -E = 2J(1- cos(4ttM/L)) 



(6) 



and consider the zero momentum | 0) as initial state. 
Note that in contrast to [l6, 27 1 we add the parameters 
M and <fi to the driving potential in Eq.©. Parameter 
M < fcmax/2 allows for the coupling of the initial state 
to very high momentum states | M) and | 2M) , and the 
parameter <f>, key to our model, allows for the coupling 
between the | s&) and | ak) basis states. 

Our choice of a two-harmonic spatial potential times 
a monochromatic time-dependent potential implies that 
symmetry St (labelled 5*2 in 0) is not broken for t s = 
tt/(2u). Therefore (p)jj = and no asymptotic current 
is possible for our system. The average current generated 
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FIG. 2: (color online) Sketch of the second order process 
in Eq. (|7J that takes place for weak coupling V/J < 1 be- 
tween the three resonant Floquet states. Dashed horizontal 
line shows the resonance quasienergy for zero driving strength 
and 5 indicates the quasienergy shifts induced by the poten- 
tial in Eq. (0 ■ Arrows represent the couplings between the 
Floquet basis states induced by the driving. The optimal cur- 
rent shown in Fig. [T]is obtained when the mixing between the 
symmetric and antisymmetric states is maximized and both 
couplings in this diagram have the same weight. 



after n cycles arises only from crossed terms between the 
cyclic states. Our aim is to maximise the average cur- 
rent during any experimental time t e . The key ingredi- 
ents are to keep few terms in the sum in Eq.Q, with 
small exponents and relevant prefactors. The two first 
are achieved by tuning the resonance in Eq. ^ with 
weak driving V/J < 1. We show that the prefactors 
can be successfully optimized if the quasiresonant cyclic 
states that have non-zero projection into the initial zero 
momentum state mix the symmetric and antisymmetric 
momentum states, which is obtained for (f> ^ for / 
integer. For cf> = lir, accidental degeneracies could in 
principle allow to obtain a small non asymptotic current 
for some specific parameters and a particular value of 
the coupling [l8[. In contrast, we show that for 4> ^ ?tt/2 
one can tailor an interference between two paths of the 
same perturbative order and find optimal parameters M, 
4> and a for any L and J set by Hq, to obtain an aver- 
age current that can be tuned up to near optimal value 
X(t e ) ~ M ~ & max /4, for a time interval [0,i e ] where 
t e can be independently tuned by adjusting the driving 
strength V/J. 

We show in fig [Iji) the average current I and the os- 
cillating average current per cycle function of 
the number of cycles. We note that the average current 
achieves a maximum M in recoil units and, as expected, 
vanishes for long times. We observe in fig QJ)) that the 
current changes direction with a sign change in (f> and t e 
scales with V~ 2 lo. For the weak driving strength used 
here the current is nearly zero (I < 10~ 3 ) for <f> = 0. 

The previous results arc obtained from a full numer- 
ical calculation. In order to further understand the ef- 



fect we follow the methodology of |27| and develop an 
approximate perturbative model. We find that for our 
potential this simplified model also explains the main 
features and gives very accurate results for low driving. 
Close to the resonance Eq. ^ and for weak driving 
V / J < 1 the dynamics of the system involve only three 
Floquet basis states {| s^m, 2) , | 0, 0) , | a,2M, 2)}, where 
\j) = \k,n) with (t\n) = e - inu)t . We apply time- 
independent perturbation theory in Floquet space, using 
the T-matrix approach T(e) = V + VGo(e)T [281, where 
^o(e) = J2j ^e-e° an< ^ £ j = Ek — nu. Around the ground 
state quasienergy e = £q, the first non-zero term connect- 
ing the three states is given by the second order in the 
expansion T(£q) — VG^e^V which reduces to 
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(7) 



A sketch of the relevant processes is depicted in Fig|2] 
and the exact values of the matrix elements, inverse of 
quasienergy differences which depend on M/L, can be 
found in the complementary material. The energy shifts 
arc 8 and the couplings D, that correspond to each part 
of the potential in Eq. (0 arc indicated by the subindcx 
M or 2M. Remarkably, O2M oc a 2 sin(2</)) and thus the 
coupling between the symmetric and antisymmetric basis 
states requires 0^0. The other effect of <j> is to bring 
those states closer in energy 8 s — 5 a oc a 2 cos(2tp). The 
optimal current is obtained when the cyclic states (re- 
lated to eigenvectors of the above matrix) mix the three 
basis states on an equal foot, corresponding to the two 
second order processes sketched in Fig. [5] being of the 
same order. Due to the structure of the spectrum, this 
optimal mixing can be reached for M < fc max /4. Then all 
the matrix elements in Eq. ([7]) are of order O ~ 1/cj and 
the quasienergies Ej ~ V 2 /ui, leading to the time scale of 
the dynamics shown in Fig. [TJd). 

For an initial state | ip(0)) = \ 0), the average current 
at cycle m after evolution with the effective Hamiltonian 
Eq. ([7]) reduces to a sum of 3 oscillatory terms with dif- 
ferent frequencies and the same weight 



l m = C(a, 



j<j' 



sm^nTAejj'), 



(8) 



where C(a, 4>) = 4ic^c 2 (p)i2 with p = 2(| s 2 m) {a 2 M \ + 
\o,2m)(s2m\) m the reduced basis and Ae^y := 
( — (ej — ej'). Once the prefactors c\c 2 are opti- 
mized for any given M < fc max /4, the current is linear 
with M. Thus, we can set M opt as the closest integer 
to fcmax/4 and obtain a robust near optimum value for 
the integrated current of X ~ M opt . As shown in figure Q] 
(dashed line), the three-mode approximation Eq. (JSj) fits 
perfectly the exact numerical results. 

We show the current amplitude C in the left panel of 
figure [3] for different parameters a and (f> . As explained 
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FIG. 3: (color online) Left panel: Amplitude C/AM of the 
current in Eq.((8]) for different values of the potential parame- 
ters cj> an d a for an initial zero momentum state and L — 41, 
M — 5, J/co — 1.0398. Right panel: Eigenenergies of the 
effective T-matrix in Eq.© in units of V 2 /(2w) as a function 
of <j> for q = 1.2. For fixed values 4> an d a the average current 
per cycle in Eq. (j8j) is a sum of sinusoidals with an ampli- 
tude C shown in the left panel and frequencies given by the 
quasienergy differences. 



above, it attains its maximum at a op t — 1.2 when both 
terms of the driving potential have the same weight , it is 
periodic in <p with 7r periodicity and has vanishing values 
for <f> = lir/2 with I integer. For (j> opt ~ ±7r/4 one can see 
in the right panel of fig. [3] that the quasienergies become 
equidistant and thus there is only one relevant energy 
scale Aef = 0.205V 2 /u) in the sum in Eq. ([5]) while the 
other sinusoidals oscillate with half this frequency. One 
can then average I m over different periods to obtain T 
which achieves its maximum Z(i c = n c T) ~ M after 
n c = 1.81 (V/u>)~ 2 cycles as shown in fig. [TJb). 

In the context of cold atoms it may be of interest not 
only the generation of a current from an initial zero mo- 
mentum state, but also the control of the quantum state 
of the system. We plot in fig. H] the particle state in 
the momentum basis and the average current per cycle 
and the average kinetic energy f Q T dt(H )/T. We observe 
that the zero momentum state can be indeed converted 
into an almost pure momentum state | ±2M). One could 
then switch off the driving, thus breaking the time re- 
versal symmetry St, and use this scheme to generate an 
asymptotic current. This is an example of the high con- 
trollability of our system. 

Finally, let us analyze the feasibility of the model. 
We can summarise our findings in a simple recipe. Hq 
sets the energy scale J and the length L of the sys- 
tem. We can obtain a final average current I of nearly 
(L — l)/8 in recoil units for a given time interval [0,t c ] 
if one tunes a driving potential in Eq. (|5|) with param- 



eters us from Eq. 



Mop 



pt; ^opt ) ^opt 



defined above 

and V = ^11.377(1 - cos(47rM op t/£))/i c with the con- 
straint that V I J < 1. We show in fig[3]that small changes 
in a and <f> around optimal values only slightly affect the 
current. Smaller M would reduce the maximum average 
current attained and require adjustment of the resonance 
condition in Eq.©. Thus the only actual requirements 
of our model are that the system is tuned to resonance 
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FIG. 4: Upper panel: Real (thick bar) and imaginary part 
(thin bar) of the {(-2M | , (0| , (2M |} | ip(t)) at t = and at 
the times showed by vertical lines in the lower panel. Lower 
panel: Average current per cycle T m /M in recoil units and 
average tunneling energy per cycle (dashed line) in units of 
1/uj as a function of time. Same parameters as in fig. [JJa). 
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FIG. 5: Average current in recoil units in Eq. ((3| as a function 
of time for an initial zero momentum state. Parameters M = 
5 , L = 41, J = 1.0398, </>/tt = 0.2579 and a = 1.2. Dashed 
line corresponds to V = 0.2 and ui — 1.01, thick line V — 0.5 
and u = 1, dash-dot line: V = 0.2 and u = 1 with an initial 
state (A#(0)} ~ exp(-£; 2 /L5 2 ) 



and that the driving field is weakly coupled. To show the 
robustness of the method we plot in figure O the average 
current when the system is not perfectly tuned to res- 
onance and when the driving field is stronger. We find 
that couplings up to V/ui = 0.5 and errors of 1% in the 
resonant frequency still give rise to high particle currents. 
Note that due to the resonance condition at low driving 
the effect is highly selective. If the initial state is a nar- 
row wavepacket centered at k ~ only this component 
is mixed to k = 2M, whereas other k < M components 
remain uncoupled. The reduced average current will be 
just proportional to the weight of k = state, see fig. [5] 
Wc have presented a model system where one can ob- 
tain currents with amplitudes orders of magnitude larger 
than those observed in recent experiments with coher- 
ent ratchet currents [HI]. The oscillation period of the 
current can be controlled by the amplitude of the driv- 
ing potential and in particular, by decreasing the driving 
strength one obtains currents which do not decay during 
the lifetime of the experiment. This effect is obtained 
with a potential which does not break the time-reversal 
symmetry and is due to crossed terms between the cyclic 
states of the system. The proposed scheme requires that 
the system is tuned to resonance and that the driving 
potential is weakly coupled such that only a few cyclic 



5 



states are involved in the dynamics. We have checked 
the robustness of the method by tuning the system out 
of resonance and increasing the coupling strength. Fur- 
thermore, we have shown that it is possible to control 
the quantum state and the amount of kinetic energy in 
the system, using the proposed scheme to convert a zero 
momentum state into a state with high finite momentum, 
and viceversa. 
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